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Abstract 

We analyze the recent NA48 data for the reaction Kl tt^jj with and without 
the assumption of vector meson dominance (VMD). We find that the data are well de- 
scribed by a three-parameter expression inspired by 0{p^) chiral perturbation theory. 
We also find that it is impossible to fit the shape of the decay distribution and the 
overall rate simultaneously if one imposes the VMD constraints on the three parame- 
ters. We comment on the different fits and their implications for the CP-conserving 
component of the decay Kl 7r^e'^e~ . 



1 Introduction 



In Ref. we examined the KTeV data for the mode Kl — > 7r°77 using a more general 
(three-parameter) description than the one used by KTeV. The latter has become the norm 
in the literature and it follows from an assumption of vector meson dominance (VMD) ||^ in 
conjunction with the parametrization inspired by 0{p^) chiral perturbation theory of Ref. 
We argued that VMD in this decay is an experimental question and, therefore, that it should 
not be an input to the data analysis. We found a least squares best fit to the data within 
our approach that was slightly better than the usual fit. However, it was hard to reach 
definitive conclusions because the necessary information is not made available by KTeV. 
Nevertheless, we motivated our more general approach by showing that there are important 
contributions to this decay from intermediate /2(1270) mesons that do not conform to the 
VMD parametrization 

In this paper we present our three-parameter fit for newly released data from NA48 
Q. This is important for the following reasons. First, the NA48 data are significantly 
different from the KTeV data and leads to different conclusions regarding the CP-conserving 
contribution to 7r°e"*"e^ 0^0 ■ Second, NA48 has presented their data in a form that 

allows us to directly compare our general fit to the usual VMD fit. This allows us to show that 
whereas it is possible to fit the decay distribution dT/dm^^ equally well with the general and 
VMD approaches, only the former is capable of fitting simultaneously the decay distribution 
and the total decay rate. 

We also discuss two additional issues. First, we show that the two fits to the decay 
distribution (the general and the VMD parametrizations) correspond respectively to con- 
structive and destructive interference between two amplitudes. Second, we comment on the 
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dependence of the fit on the parameter Og, which is extracted from K tttttt decays and 
which has a large uncertainty (jTUl ITT . 



2 Parametrization of the data 

The Kl — > '7r°77 amphtude in the hmit of CP violation can be written in terms of two 
independent invariant amplitudes, A and B M, 
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The Fermi constant and the Cabibbo angle are included in the overall constant Gg = 9.1 x 
10^^ GeV~^ and a^^j ~ 1/137 is the usual electromagnetic fine structure constant. To 
parametrize these amplitudes in a form inspired by Oijf') chiral perturbation theory and 
dispersion relations Ref. H| proposed the use of: 
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where the dimensionless kinematic variables are 
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and the scale of chiral symmetry breaking is Airf^ ^1.17 GeV. 

The form in Eq. (H) does not correspond to a complete calculation at order in 
chiral perturbation theory. Rather it contains the complete one-loop calculation at order 
|jl2[ and two additional ingredients containing some corrections of order p^ |jl3|, The 
non-analytic terms in Eq. that multiply the factors 05 and a^(z) attempt to incorporate 
the strong rescattering in the two-pion intermediate state that occurs at one-loop. They 
arise from the inclusion of corrections to the K ^ Sti amplitudes [l^. The values 
of and 03 are extracted from data and the functions F{z) and R{z) can be found in the 
literature [Q. The three constants a^, ctg and f3 originate in counterterms appearing in the 
p^ weak chiral Lagrangian 0]. 



The analysis of ^ 37r in Ref. |10| indicates that 



a^{z) = 0.38 + 0.13 Yq - 0.0059 , = 6.5 , 



Yn 



K 
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A very recent analysis of if 37r data results in |TT 



02 = 6.8 ±2.4 . 
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In the analysis of Ref . ^ , which has become standard, the three unknown constants 
were fixed in terms of the contribution they receive from vector-meson exchange, supple- 
mented with a minimal subtraction ansatz: 

Oil = -^ay , 

a2 = 12ay — 0.65 , 

P = -Say -0.13, (6) 

and this form has been used both by KTeV and by NA48 to fit their data. In Ref. [Q] we 
argued that this ansatz imposes a correlation on f3 that is not desirable for a prediction of the 
CP- conserving contribution to —>■ Tx^e^e~ . With the new NA48 data we can go further 
and conclude that the VMD ansatz does not provide a good description of 7r°77. 

3 Fitting the shape of the dT/dm^^ distribution 

NA48 has recently released their result for 7r°77 [0. They chose to analyze their 

data using Eq. (|^) with the VMD assumption, and they found ay = —0.46. To obtain 
this number they fit the shape of the distribution dV / dydm^^ without attempting to fit 
the branching ratio. NA48 has published in Table 2 of Ref. |^ sufficient information to 
fit the distribution dV/dm^^. They present the number of unambiguous events, estimated 
background and acceptance for each 20 MeV bin in m^^. 

We begin our analysis with a fit to the shape of the dV/dm^^ distribution, ignoring 
the measured branching ratio, to compare with the fit performed by NA48. We do this both 
using the VMD assumption and with the general approach. We calculate the number of 
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events predicted in each bin as 



Ni = N 



\(Jjl / L-yy / 
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Acceptance, + Background^ , (7) 



where is a normahzation chosen to match the total number of events and N{Kl) = 23.9 
X 10^ is the number of decays in the fiducial volume. The arbitrary normalization allows us 
to fit the shape of the distribution while ignoring the overall rate. 

We use data from 17 out of 20 bins presented in Table 2 of Ref. 0. We exclude two 
bins in the m^^ region near the vr" mass which do not have any events due to kinematic cuts, 
and we also exclude the last bin with no events because it lies outside the physical region. 
We perform a least squares fit using Poisson statistics for the bins with small number of 
events following Ref. WE . 



With this procedure, and the VMD ansatz, we reproduce approximately the NA48 
best fit. We obtain ay = —0.455 with a y^/dof = 18.5/16. We show this result in Fig. |l] 
where we superimpose our best three-parameter fit which has a x^/dof = 14.6/14. The two 
fits are nearly identical as can be seen in the figure and they are indistinguishable statistically. 
Nevertheless, when they are both expressed in terms of the three general parameters one can 
see they correspond to very different solutions. For the general fit, 

= 4.57 , = -3.89 , p = 0.75 ; (8) 

whereas for the VMD fit (in terms of a^), 

= 1.82 , ^2 = -6.10 , p = 3.51 . (9) 

For the case of the three-parameter fit we find that a-^ and 02 correlated as was 
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Figure 1: Two different fits to the data from Ref. 0, as explained in tlie text. Tlie solid 
line is a one-parameter fit corresponding to Eq. (|^), the dashed line is the three-parameter 
fit shown in Eq. (H). 

discussed in Ref. [Q, so that there are many other fits with a near the minimum for the 
same value of /3. 

As stated above, neither one of these fits reproduces the experimental rate, B{Kl 
^°77) = (1.36 ± 0.03 ± 0.03) x 10^^ 0. The theoretical branching ratio predicted for 
ay = —0.455 (the NA48 value) is B{Kl — > 71^77) = 1.1 x 10"^, and the one predicted for 
the three parameters in Eq. (H) is B{Kl TT^'-f-f) = 1.0 x 10~^. 
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It is instructive to show the three separate contributions that result from Eq. (|^) 
to the differential decay rate dB{KL Tt^'-yj) / dm^^. The three terms correspond to the 
absolute square of the A and B amplitudes and to their interference, |y4p, and Re(y4*i?), 
respectively. We show these quantities in Fig. ^ for the best three-parameter fit and Fig. ^ 



for the best ay fit. 
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Figure 2: The contributions from \A\'^ (solid line), (dashed line), and Iie{A*B) (dot- 
dashed line) are plotted vs. the invariant two-photon mass m^^ in terms of the number of 
events for the best three-parameter fit. 



In both of these figures the solid line represents the contribution from the dashed 
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Figure 3: Same as in Fig. ^for the best ay fit. 

line the contribution from \B\'^ and the dot-dashed fine the interference. We observe that the 
three-parameter fit corresponds to constructive interference between the A and B amphtudes, 
whereas the Oy fit corresponds to destructive interference. Unfortunately it appears that it 
is not possible to determine experimentally the sign of this interference. However, as we 
show below, the total rate for the process discriminates between the VMD ansatz and the 
general form of the amplitude. 
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4 Simultaneous fit to the shape of the dT/drrij^ distri- 
bution and to the decay rate 



To obtain a fit that reproduces the observed branching ratio we proceed as in Eq. (|^) but 
removing the arbitrary normahzation, 



Acceptance^ + Background^ , (10) 



with the same notation of Eq. (0). We first attempt this fit with the VMD ansatz and 
find that it is impossible to obtain a good fit. Our least squares fit using the VMD ansatz 
occurs for ay = —0.63 and has a x^/dof = 74.8/16. We show this result as the solid line in 
Figs. ^ and ^ The implied branching ratio is B{Kl — > 7r°77) = 1.25 x 10^^ and ay = —0.63 
corresponds to 

ai = 2.51, a2 = -8.19, /3 = 4.89 . (11) 

Our best three-parameter fit, on the other hand, has a /dof = 18.5/14 and is shown 
as the dashed line in Figs. §and ^ It implies a branching ratio B[Kl 71^77) = 1.36 x 10~^ 
in good agreement with the measured one. The parameters for this best fit are, 

= -2.46 , ^2 = -2-51 , (3 = 0.16 . (12) 

We conclude from Fig. ^ that the VMD ansatz cannot reproduce the shape of the spec- 
trum and the total decay rate simultaneously, but that the general formula, Eq. (0) does 
accommodate both. 

For completeness we show in Figs. ^ and the theoretical dT/dy distributions for both 
the tty result from Eq. (|TTp and the three parameters given in Eq. (|r^). Fig. |^ is restricted 
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Figure 4: A simultaneous fit to the shape of dT/dm^^ and to the decay rate. The solid line 
is a one-parameter fit corresponding to Eq. (|ll|), the dashed line is the three-parameter fit 
shown in Eq. (|12D. 



to events with m^^ < 0.24 GeV/c^. There are no data available in this form, so at this point 
we are not able to perform a fit and we can only present our predictions. We point out that 
the three-parameter fit yields a flatter distribution than the ay fit. 
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Figure 5: An enlargement of Fig. ^ for m^^ < 0.24 GeV/c^. 
4.1 Dependence on a2 

We now consider the dependence of our results on the parameter 03 that appears in the B 
amplitude. This parameter is extracted from K decays and up to now we have used 

the value 03 = 6.5 [Q. However, the value of this parameter has a large uncertainty, of order 
~ 35 %. For example, from the recent analysis of Ref. ||Tl|| one extracts 03 = 6.8 ± 2.4. 



The analytic form for the B amplitude in Eq. clearly indicates that and (3 
are correlated and this is confirmed by our numerical study. It is possible to obtain many 
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Figure 6: Theoretical dV/dy distributions. The sohd and dashed hues are predicted using 
as input the results given in Eq. ( pA]) and Eq. (0), respectively. 

equally good fits to the data with different values of and (3. For example if we take the 



central value from Ref. [ITT] and 1-sigma deviations from it, we find good fits to the shape and 
spectrum with the values listed in Table |I[ This is not possible with the ay parametrization, 
where we cannot find a good fit for any of these values of 
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Figure 7: Same as Fig. | but with a cut at m^^ < 0.24 GeV/c^. 

5 CP-conserving contribution to Kl — > Tr^e^e" 

We now turn to the estimate of the CP-conserving contribution to vr^e'^e" using 

the model of Ref. [Q. Using the results of the fit to the shape of the distribution only, 
Eqs. (g) and (g), we find 



4.0 X 10 vector meson dominance 
2.0 X 10^^'^ three-parameter fit . 
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Table 1: Three-parameter best fits for tliree different values of 02, corresponding to its central 
value from Ref. [|rT| and its 1-sigma deviations. 

Notice that these two numbers are an order of magnitude smaller than what is obtained using 
the KTeV data instead (see Eq. (11) of Ref. ]1|1)- We can see from Fig. || why the NA48 
result implies a much smaller Bcpc{Kl vr^e+e^) than the KTeV result (/? = —5 
for the three-parameter fit or /? = 7.5 for the ay fit). These two points are shown as the two 
internal dotted lines in Fig. p. It is clear from this figure that the NA48 results correspond 
to a Kl — ^ TT^'-ff that produces a minimal CP-conserving contribution in Kl —* vr^e^e", i.e. 
it indicates that the two photons have a negligible D-wave component. The VMD result in 
Eq. (|T3p is consistent with the result reported by NA48. The latter is based on an analysis of 
the low m^^ region only and yields Bcpc{Kl vr°e+e-) = (4.71?:^) x 10"^^ 1- The NA48 
result is obtained from data with m^^ below 110 MeV and is therefore model independent 
because in that region the B amplitude dominates and the correlation with the A amplitude 
implied by the VMD ansatz disappears. 

If we use the results of the fits to both rate and spectrum, Eqs. ([TTp and (0), we find 
instead, 

jj , + _ J (13.8l2i) X lO^^'^ vector meson dominance ^ 

BcPciK^ - vr e e ) - I ^^ ^+,.7^ ^ ^^^,3 ^^ree-parameter fit . 

These two points are shown as the external dotted lines in Fig. ^ Not surprisingly, the 
general three-parameter fit continues to agree with the model independent NA48 limit as it 
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Figure 8: CP-conserving contribution to vr^e^e" as a function of (3 with 02 = 6.5 \V\. 

The dashed hne shows the absorptive contribution and the sohd hne the model of Ref. 
The enlargement shows the results for the branching ratio vs. the four values of /3 = 0.16, 
0.75, 3.51 and 4.89 from the three- and one-parameter fits discussed in the text. These are 
marked by vertical dotted lines. 

gives a good fit to both the rate and spectrum. On the other hand, the fit in terms of ay 
alone does not reproduce the data very well and we can dismiss its implication of a larger 

In Fig. § we see why there are two different solutions for (3 that result in the same 
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Bcpc{Kl 7r°e+e~). This CP-conserving component depends quadratically on the B{z) 
amphtude of — > 7r°77, and therefore there are two values of [3 for any given Bcpc{.Kl — > 
TT^e^e~). They correspond to constructive and destructive interference between the term 
with 02 and /3 in Eq. (|]). 

6 Conclusion 

We have shown that the NA48 data for the reaction Tc^'j'y can be accommodated 

nicely by the theoretical expression based on chiral perturbation theory at order p^. With 
this expression it is possible to describe simultaneously the total rate and the shape of the 
spectrum, which is not possible with chiral perturbation theory at order p'^ We have 

also shown that the commonly used VMD ansatz fails in this case, and that it is impossible 
to fit both the rate and the shape of the spectrum if this ansatz is adopted. 

We have also shown that it is possible to obtain a good fit to this mode for different 
values of the poorly known parameter from K ^ 3tt decays. This indicates both that 
Kl — »• 7r°77 cannot provide additional information on the value of 03, and that not knowing 
its precise value does not affect our ability to describe the features of Kl 7r°77. 

Although we do not have sufficient information to perform a similar comparison for 
the KTeV data, we note that the value of ay reported by KTeV 0, ay = — 0.72±0.05±0.06 
predicts a branching ratio B{Kl Tr^'j'-f) = (1.36 ± 0.06 ± 0.07) x 10^^ in conflict with the 
measured value, B{Kl vr077) = (1.68 ± 0.07 ± 0.08) x 10"^. 

The new results from NA48 indicate a very small D-wave component for the photon 
pair and this leads to a prediction of a negligible CP-conserving background to Kl 
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7r^e~^e . We have shown that this result is not an artifact of the VMD ansatz and that it 
holds in the general parametrization. This result is at odds with the earlier KTeV data and 
we must wait for the new KTeV results to see how this discrepancy is resolved. 



Acknowledgments 

This work was supported in part by DOE under Contract Number DE-FG02-01ER41155. 
We thank R. Wanke and M. Martini for helpful discussions of the NA48 results. We thank 
James Cochran for bringing Ref. [1T3I to our attention. 



References 

[1] F. Gabbiani and G. Valencia, Phys. Rev. D 64, 094008 (2001) ||hep-ph/010500l . 
[2] KTeV Collaboration, A. Alavi-Harati et al. Phys. Rev. Lett. 83, 917 (1999) [flep^ 



ex/9902029|| 



[3] G. Ecker, A. Pich, and E. de Rafael, Phys. Lett. B 237, 481 (1990); L. M. Sehgal, Phys. 
Rev. D 41, 161 (1990); P. Ko and J. L. Rosner, ihid. 40, 3775 (1989); P. Ko, ihid. 41, 
1531 (1990). 

[4] A. G. Cohen, G. Ecker, and A. Pich, Phys. Lett. B 304, 347 (1993). 

[5] Other contributions are also discussed in J. Bijnens, S. Dawson, and G. Valencia, Phys. 
Rev. D 44, 3555 (1991); S. Fajfer, ihid. 51, 1101 (1995) ||hep-ph/94063T7|] ; Nuovo Ci- 
mento A 109, 1439 (1996). 



[6] NA48 Collaboration, A. Lai et al. |hep-ex/0205CTTU 



17 



[7] J. F. Donoghue, B. R. Holstein, and G. Valencia, Phys. Rev. D 35, 2769 (1987); L. M. Se- 
hgal, tbid. 38, 808 (1988); R Heiliger and L. M. Sehgal, ibzd. 47, 4920 (1993); J. Kambor 
and B. R. Holstein, iMd. 49, 2346 (1994) ||hep-ph/93l0^ ; Some reviews of this and 
other kaon decay modes are J. L. Ritchie and S. G. Wojcicki, Rev. Mod. Phys. 65, 
1149 (1993); L. Littenberg and G. Valencia, Annu. Rev. Nucl. Part. Sci. 43, 729 (1993) 



hep-ph/ 93032251 ; G. D'Ambrosio and G. Isidori, Int. J. Mod. Phys. A 13, 1 (1998) 



|hep-ph/961128l ; A. R. Barker and S. H. Kettell, Annu. Rev. Nucl. Part. Sci. 50, 249 



(2000) [|hep-ex/ 0009024 



[8] G. Ecker, A. Pich, and E. de Rafael, Nucl. Phys. B303, 665 (1988) 



[9 
[10 

[11 
[12 

[13 
[14 
[15 
[16 



J. F. Donoghue and F. Gabbiani, Phys. Rev. D 51, 2187 (1995) ||hep-ph/9408'Ml 



J. Kambor, J. Missimer, and D. Wyler, Phys. Lett. B 261, 496 (1991). 
J. Bijnens, P. Dhonte, and F. Persson, piep-ph/0205341. 



G. Ecker, A. Pich, and E. de Rafael, Phys. Lett. B 189, 363 (1987); L. Cappiello and 
G. D'Ambrosio, Nuovo Cimento A 99, 153 (1988). 

L. Gappiello, G. D'Ambrosio, and M. Miragliuolo, Phys. Lett. B 298, 423 (1993). 
J. F. Donoghue, E. Golowich, and B. R. Holstein, Phys. Rev. D 30, 587 (1984). 
G. J. Feldman and R. D. Cousins, Phys. Rev. D 57, 3873 (1998) ||physics/971102T 



NA31 Collaboration, G. D. Barr et al. Phys. Lett. B 242, 523 (1990); 284, 440 (1992); 
328, 528 (1994); E731 Collaboration, V. Papadimitriou et al Phys. Rev. D 44, 573 
('1991). 



18 



